Confinement-induced Resonances in Quasi-one-dimensional Traps with Transverse 

Anisotropy 
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We study atom-atom scattering in quasi-one-dimensional geometries witli transverse anisotropy. 
By assuming an s-wave pseudo-potential of contact interaction, we show that the system would 
exhibit a single confinement-induced resonance, where the scattering process degenerates to a total 
reflection as a one-dimensional gas of impenetrable bosons. For a general form of interaction, we 
present a formal calculation based on the two-channel model and conclude the existence of only one 
confinement-induced resonance. Our findings are inconsistent with a recent experiment by Haller et 
al. [Phys. Rev. Lett. 104, 153203 (2010)], where a splitting of confinement-induced resonances has 
been observed in an anisotropic quasi-one-dimensional quantum gas of Cs atoms. 
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I. INTRODUCTION 

The interest on low-dimensional systems has been re- 
cently intensified by the experimental development of 
manipulating ultracold quantum gases in optical lat- 
tices [ll-Q and on atom chips . With the aid of tun- 
ing an external magnetic field through a Feshbach reso- 
nance, these techniques provide a fascinating possibility 
of investigating distinct physical properties in quasi-low- 
dimensionality with a controllable particle-particle inter- 
action. Up to now, the Berezinskii-Kosterlitz-Thouless 
(BKT) transition has been observed in a series of 
quasi-two-dimensional (quasi-2D) pancake shaped Bose 
gases [1], and the strongly interacting Tonks-Girardcau 
(TG) Bose gas has been achieved in a setup of 2D op- 
tical lattices 0, H, [l^, [HI • Of particular interest is the 
quasi-one-dimensional (quasi-lD) regime, where only the 
ground state of transverse motion is significantly popu- 
lated. This regime is important partly due to its the- 
oretical simplicity. In fact, the homogeneous system of 
ID Bose gas with ^-contact interaction is known as one 
of the few fully integrable quantum systems [l^ . In a 
finite system with infinitely strong 5 interactions, a TV- 
body Bose system has been proved to correspond via a 
one-to-one mapping with the highly correlated states of 
the corresponding noninteracting Fermi gas [l3l - [T8j . 

The quasi-lD geometry can be realized by arranging 
a 2D optical lattice in the transverse {x-y) plane, such 
that the kinetic and the interaction energy of the parti- 
cles are insufficient to transfer particles to transversally 
excited energy levels. However, even in this ultimate 
limit, one should bear in mind that the transverse de- 
grees of freedom are of great importance and could man- 
ifest themselves in many aspects of the system properties. 
For example, when the particle-particle interaction can 
support a bound state, the binding energy will serve as 
an additional energy scale, which can be comparable or 
even exceed the transverse confinement. In this situation, 
the transversally excited states will be inevitably popu- 
lated, and have to be taken into account for a proper 



description of the system [l9|, |20|. On the other hand, 
when the interaction does not support a boimd state, 
the level structure of transversally energy states will pro- 
vide extra channels for the ID scattering process. The 
presence of these channels gives a scheme to tune through 
the ID scattering resonance by either changing the three- 
dimensional (3D) s-wave scattering length, or varying the 
transverse confinement, which is known as confinement- 
induced resonance (CIR) as proposed by Olshanii [2]| . 
Around the resonance point of CIR, the effective ID cou- 
pling constant can be varied from — oo to -|-oo, and the 
ID scattering process degenerates to a total reflection, 
thereby creating a gas of impenetrable bosons. 

The underlying physics of CIR is in the same spirit of a 
Feshbach resonance between open and closed scattering 
channels (2^ . In the presence of transverse confinement, 
the scattering of atoms in the transversally ground state 
assumes the open channel, while the transversally ex- 
cited states as a whole can support a bound state and 
serve as the closed channel. According to the Feshbach 
scheme, the resonance takes place when the bound state 
of the closed channel degenerates with the continuum 
threshold of the open channel, which predicts the loca- 
tion where GIR occurs. By assuming a 3D Fermi-Huang 
pcsudo-potcntial for the particle interaction, it is shown 
that the GIR takes place in a cylindrically symmetric 
quasi-lD harmonic trap at a specific ratio of as/a± ~ 
— l/C(l/2), where Os is the 3D s-wave scattering length, 
aj_ = \/fi/ is the length scale of transverse harmonic 
oscillator with trapping frequency uj± = cux = ojy, fJ. is 
the reduced mass of the two colliding particles, and ({x) 
is the Ricmann Zeta function [2l|. Further theoretical 
investigation extends the discussion of CIR to pancake 
shaped quasi-2D geometry 



23| 



to general forms of in- 
teractions j24| - |26l |. to the three-body [l^] and the four- 
body [i^ scattering processes, as well as to a pure p-wave 
scattering of fermions [2^ . 

Experimental evidence of GIRs has been recently re- 
ported for bosonic 0, [ll|, |33| and fermionic @ sys- 
tems. In particular, an experiment by Haller et al. per- 
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FIG. 1: Locations of CIRs as a function of the transverse 
anisotropy rj = ujx/ijJy By assuming a 3D Fermi-Huang 
pseudo-potential, we find a single CIR for all values of trans- 
verse anisotropy (solid). For small anisotropy of ry « 1, the 
resonance position shifts downwards with increasing at a 
slope of —1/4 (dashed). 



that the sphtting of peaks of atom loss and heating rate 
observed in the experiment cannot be simply understood 
by the theory of CIR. 

The remaining of this manuscript is organized as fol- 
lows. In Sec. ini we discuss the s-wave scattering process 
of two atoms in a quasi- ID geometry with transverse 
anisotropy. By assuming a (S-contact pseudo-potential, 
we calculate the position of CIR for a wide range of trans- 
verse anisotropy. This result is then confirmed by another 
calculation of the closed channel bound state energy, as 
discussed in Sec. IIIII Considering the fact that the CIR 
takes place at the exact position where the bound state 
of the closed channel degenerates to the threshold of the 
open channel, we can obtain resonance condition from a 
different aspect. In Sec. IIVI we extend the discussion to 
a general form of interaction, and present a formal the- 
ory for two-atom scattering within a two-channel model. 
Finally, we summarize our conclusions in Sec. |Vl We put 
some detail calculations for Sec. IIVI in Appendix[Dl while 
some necessary calculations for Appendix [D] arc posted 
in Appendix EEl 



forms a systematical investigation on CIRs in a ultra- 
cold quantum gas of Cs atoms confined in a quasi- ID ge- 
ometry with transverse anisotropy [soj . By varying the 
3D s-wave scattering length around a magnetic Feshbach 
resonance, they determined the positions of CIR from 
atom loss and heating rate measurements. Of particu- 
lar interest is the CIR was found to split into two res- 
onances upon introducing anisotropy to the transverse 
confinement Ux 7^ One of the resonances, CIRl, 

shows a pronounced shift to higher values of as /ay by 
increasing anisotropy rj = Ux/ujy. The second resonance, 
CIR2, shifts instead to lower values of as/ ay. Here, 
ay = y/h/J^jLLUy) is the harmonic oscillator length in the 
y-direction. 

In this manuscript, we present a theoretical investiga- 
tion on the CIR in a quasi-lD harmonic trap with trans- 
verse anisotropy. By using a 3D Fermi-Huang pseudo- 
potential, we show that only one single CIR would be 
present for all values of anisotropy ratio rj. For small 
transverse anisotropy, the resonance position of as/ ay 
shifts downwards with increasing rj. This is in the same 
trend as the experimental observation of CIR2 [s^ , but is 
systematically above the measured values. The appear- 
ance of a single CIR can be well understood by notic- 
ing that the s-wave contact interaction can only support 
one bound state in the closed channel of excited trans- 
verse modes, and the resonance can only occur when this 
bound state energy coincide with the threshold in the 
open channel of ground transverse mode. In order to 
further validate the observation to a general form of s- 
wave interaction, we imply a formal scattering theory of 
a two-channel model, and draw the same qualitative con- 
clusion of a single CIR. Upon these findings, we conclude 



II. HAMILTONIAN AND TWO-BODY 
SCATTERING PROBLEM 



We consider the problem of two interacting atoms with 
mass m colliding in a 2D harmonic confinement in the 
transverse x-y plane 
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(1) 



In general, the transverse trapping frequencies uix,y along 
the X and y directions are different, with transverse 
anisotropy ratio rj = tOx/ojy. The atomic motion along 
the z-axis is assumed to be free. Then the Hamiltonian 
of the system is given by 

H = -^Vl--^Vl+Uir,) + Uir2)+Vir,-r2), (2) 
Zm 2m 

where ri and r2 denote the positions of the two atoms, 
respectively. V{ri — is the interaction potential be- 
tween the two atoms. For sufficiently low energies, the 
3D atom-atom scattering is dominated by s-wave scat- 
tering. Therefore, in the present and the next sections, 
we model the interactionpotential V hy a. Fermi-Huang- 
Yang pseudo-potential |3l| in the form of regularized 5- 
function 



ATTfi^as d 

Viy) = S{r)—r, 

m or 



(3) 



where as is the 3D s-wave scattering length. The dis- 
cussion of a general form of interaction will be left in 
Sec.IVl 

One great advantage of harmonic trapping potential is 
the center-of-mass (COM) and relative degrees of free- 
dom can be decoupled. Substituting new coordinates 
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r = Ti — r2 and R = (ri + V2)/'2., we separate the full 
Hamiltonian into the COM part Hqom and the relative 

part ifrel 



-f^CO 



M 



— V^ + 2C/(R), 



Hr, 



2M 

^ 

2/i 



C/(r) 



y(r), 



(4) 
(5) 



where M = 2m and fi = m/2 are the total and re- 
duced masses, respectively. The eigenfuctions of the 
COM Hamiltonian is simply a combination of harmonic 
oscillators in the transverse plane and plane waves along 
the z-direction, which can be safely dropped out from our 
discussion. The relative Hamiltonian can be expressed in 
a dimensionless form, in which all lengths are scaled in 
units of Uy = y/h/Jpujy), and all energies are scaled in 
units of fkdy. The eigcnfunctions 4'(r) of the relative mo- 
tion thus are governed by the Schrodinger equation 



1. 
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d 



Vr + ^iv + r) + 2TTa,5{r)—r 



*(r) = £;*(r), 
(6) 



where the eigenenergy E = (/c^ + ^7 + l)/2 is the summa- 
tion of the longitudinal kinetic energy fc^ and the ground- 
state energy of the transverse confinement. In this paper 
we assume the incident wave function is in the transverse 
ground state and longitudinal kinetic energy is smaller 
than energy spacing between the ground and first ex- 
cited transverse mode, i.e., fc^ < minjTy, 1}. In this case 
the asymptotic scattering-state wave function in the limit 
|z| — >■ oo would be definitely in the transverse ground 
state. 

As a result of the zero-range pscudo potential ^(r) in 
([3]), the scattering state with odd parity along the z di- 
rection does not get any phase shift. Then the relevant 
scattering amplitude is zero. Thus, we we focus on the 
scattering state with even parity in the following discus- 
sion. 

In the region r 7^ 0, the pseudo-potential ^(r) takes 
zero value. Therefore the atomic motion with respect 
to different cigenstates of the transverse confinement are 
decoupled with each other. Then the eigen-equation © 
can be solved analytically. The solution with even parity 
and out-going boundary condition in the z direction takes 
the form 

^ix,y,z) = (cos(fc,z) + /ee*'=^l-l)7//4<i>o(V^x)$o(2/) 



where Hn^,n = \/'2nx'>] + 2ny — k^. and the summa- 
tion runs over all combinations of {nx,ny) excepts the 
{n.j; ~ 0,ny = 0) ground mode. In the expression above 
ri^^^^n{\/vx) and $n(y) are the normalized wave func- 
tions of the n-th excited states of the confinements in the 



X and y direction, respectively, with the dimension- 

less wave function of harmonic oscillator 



Hn{t). 



(8) 



Here, H„{t) is the Hermite polynomial. 

The scattering amplitude fe and the parameters Bn^uy 
can be formally extracted from the connection condition 
given by the s-wave pseudo potential ^(r) at the origin 
r = 0. To this end we substitute Eq. © into Eq. ©, 
and do the integration 



lim / dz 



+ 00 



dxdyW^''S>*„Jy^x)KAy)\ (9) 



on both sides of the equation. Then we get 
27ra 



fe 



ik 
2'Ka 



r?i/4|$o(0)p5, 

77l/4$„^(0)<i>„JO)5, 



where the factor S is defined as 



S = 



•5 T/ ^' 
r 

or 



^z-^{Q,0,z) 
oz 



(10) 
(11) 

(12) 



To derive the finial expression of the scattering ampli- 
tude, we substitute Eqs. ([TU)) and pT|) into Eq. ([7]), and 
express ^(0,0,2:) in terms of the parameter S. Substi- 
tuting this result into Eq. ([T2|) . we get the equation for 
S: 



5=1 + 



27ra< 



771/4|<I,q(0)|2 5] ^1/4 1^^(0)1 



-asri^/^P{k,,i-l)S 



(13) 



where l3{kz,'i'i) = [dzzA{z,kz,r/)]^^Q+ is the regular part 
of the series 

(14) 

Solving S from Eq. (|13p . we finally obtain the quasi- ID 
scattering amplitude 



fe 



l-ik,[ai'fj-^/^+Pik,,r])]/2' 



(15) 



In the low energy limit of — > 0, the result of scattering 
amplitude /e (fT5|) can be approximated by a scattering 
amplitude 



f!ih) = 



-1 



1 + ikzCLii 



(16) 



for a ID ^-potential Vid = 5id<^(^) of coupling strength 
giD = ^fliDi where the ID scattering length takes the 
form 



aiD 



1 



(17) 



4 



The CIR takes place at the position of aio approaches 
zero, hence indicating a complete reflection of the scat- 
tering amplitude. By introducing back the physical unit, 
we obtain the CIR condition 



1 



ay ^/?7/3(0, v) 



(18) 



From this result, we can see clearly that only one single 
CIR would be present for arbitrary values of transverse 
anisotropy ratio rj. 

We stress that the partial derivative d/dz and the sum- 
mation over transversally excited states {n^^ny)' in A 
cannot be interchanged because the series does not con- 
verge uniformly, and the summation must be done first. 
This is related to the divergent behavior of the Green's 
function at small r, which is regularized by the operator 
dr{r-). Thus, to determine the positon of CIR, we gener- 
ally need to numerically evaluate the series A(z, 0, 77) for 
a set of z, and extract the coefficient /3(0, 77) through the 
fit 

A(z, 0, r,) = + m ri) + 7(0, ti)z + --- . (19) 
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FIG. 2: Locations of CIR as a function of tlie transverse 
anisotropy rj = cj^/cjy. By varying rj from the cylindrically 
symmetric case rj = 1, the CIR position shows non-monotonic 
behavior and crosses over the magnetic Feshbach resonance to 
reach the BCS regime with negative scattering lengtli (verti- 
cal lines). 



In order to do this, we rewrite the series as 



A{z,0,rj)= J2 P^^P^y , — ' (20) 



•sJSxrj + Sy 



where and Sy are non-negative integers and the sum- 
mation runs over all possible combinations except the 
ground state (s^ = 0, Sj, = 0). The coefficient is 



(2si)! 



r(s, + i/2) 



22-.(,S,!)2 V^(,5, + l)' 



(21) 



where r(a;) is the Euler Gamma function. 

Before discussing the general cases with arbitrary 77, 
we first focus on the cylindrically symmetric trap with 
77 = 1, where analytical result is available. In this case, 
we define a new summation variable s = Sx~\- Sy and get 



A(z,0,77=l) = ^. 



-2\/sz 



(22) 



To reach this expression, we have utilized the relation 
r(m-hl/2)r(77-77i+l/2) 



E 



r(r77,+ l) r(?7, -777+1) 



(23) 



for n = 1,2,3,... being positive integers. The series of 
Eq. ([22]) has been evaluated [2lj, leading to the well- 
known result of /3(0, 77 = 1) = C(l/2) Ki -1.4603 . . . , and 
the CIR condition reads 



1 



C(l/2) 



(24) 



Next, we discuss the general case of 77 7^ 1. For a trap 
with slight transverse anisotropy A7/ = 7^ — 1 <C 1, we 
can expand the function A in Eq. to linear order of 
A77. Notice that the series in uniformly convergent for 
all finite 7;, thus it is safe to interchange the summation 
and the partial derivative of 7;, leading to 



A(z,0,7/) « A(z,0,l) + A?/ 



drj 



0{Ar]^) (25) 



,7=1 



with the coefficient 
dA 



if:(-l + 2.)e-v- 



(26) 



The first term in the bracket corresponds to the same 
summation as in Eq. (|22p . and the second term can be 
evaluated as 



ize -''"^ ^-+0(z 

z 

s=l 

where we have used the relation 

00 



0+, 



E2. 

s=l 



n+ 



2 / dxe-2^ + 0(2 
Jo 



(27) 



(28) 



As a consequence, we reach the expression of A for 7^ ~ 1 
A{z 0+,0,77) 



(29) 

Substituting the expression above into Eq. p^ . we get 
the CIR condition 



1 



"C(l/2) 



1 



77 — 1 



(30) 
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which indicates that the position of CIR would shift 
downwards with a slope of — 1/4 if we introduce a slight 
transverse anisotropy by increasing r/ — uix/tOy (see Fig. 
1). This result is in the same qualitative trend as the ex- 
perimental measurement of CIR2 [soj . but with a slope 
only about half of the experimental observations. 

On further increasing or decreasing the transverse 
anisotropy, the position of CIR can be determined by nu- 
merical evaluation and fit A to Eq. (|19p . The results are 
plotted in Fig. 2 for a wide range of r]. It is shown that 
the resonance position of as/ ay acquires a non-monotonic 
behavior of rj, and diverges at a critical frequency ratio 
r]c ~ 18.5 (or equivalently ■q~^ « 0.054), where the 3D 
scattering length crosses through the magnetic Feshbach 
resonance to the BCS side of < 0. This behavior 
can be understood by noticing that by continuously im- 
posing the transverse anisotropy, the quasi- ID geometry 
will eventually crosses over to a quasi-2D system, where 
the s-wavc CIR occurs at the BCS side of the magnetic 
Feshbach resonance [1^ [25[ . 



III. TWO-CHANNEL MODEL AND BOUND 
STATE ENERGY 

In the previous section, we present a study on 
two-body s-wave scattering process in a transversely 
anisotropic quasi- ID geometry. We find a CIR as in a 
cylindrically symmetric quasi-lD trap. Reminding that 
CIR is a Feshbach resonance between the open channel 
of the transversally ground state and the closed channel 
of exited states, the uniqueness of CIR is rooted from the 
fact that only one bound state can be supported by the s- 
wave pseudo-potential. Thus, it is illuminating to study 
the bound states of the closed channel and examine the 
CIR condition from the binding energy. 

In order to extract the detailed information of the 
eigenstate structure, we add an additional harmonic trap 
along the axial z-direction with frequency w^, and take 
the limit oi lJz — > afterwards. This procedure allows 
us to take advantage of the simple form of harmonic os- 
cillators. The dimensionlcss Hamiltonian of the relative 
motion of two interacting particle takes the form 

iJrci = -^V2 + ]^{ifx^ +y^ + £2^2) + 2™,<5(r) |- (r • ) . 

(31) 

Here, we use ay = ^hj (fiujy) as the length unit and hiOy 
as the energy limit. By taking the limit of e = ujz/ujy — >■ 
0, we will reduce to the original problem with no axial 
confinement. 

Since the CIR is in fact a Feshbach resonance between 
transversally ground mode and the manifold of excited 
modes, we need to investigate the relevant bound state 
of the decoupled excited manifold and ask that when this 
bound state energy coincides with the continuum thresh- 
old of the ground transverse mode. Thus, we proceed 
by formally splitting the Hamiltonian onto "ground" (g), 



"excited" (e), and "ground-excited coupling" {ge) parts 
according to 

Hid = Hg + He + -ffge 

= PgHPg+PeHPe + iPaHPe + PeHPg),i32) 

where Pg = |00)_l(00|, A = E(„,,n„)' \nxny) j_{nxny\ are 
the corresponding projection operators, \nxny)± repre- 
sents the eigenstate of the transverse harmonic oscillators 
with quantum number Ux and ny, and the summation 
(nxjUy)' runs over all possible combinations except the 
ground mode |00)_l. 

The ground Hamiltonian Hg has a ID coordinate rep- 
resentation, which corresponds to the motion of a one- 
dimensional particle in the presence of a S barrier and 
an axial harmonic trap. In the limit of e — >■ O, the spec- 
trum of Hg is continuous for energies above the threshold 
energy EgQ = {i] + l)/2, which represents the zero-point 
energy of the transverse mode. In the same limit, the 
spectrum of the excited Hamiltonian He is also continu- 
ous for energies greater than Eeo = 3(77-1- 1)/2. But as we 
will see below, He also supports a bound state with en- 
ergy Eeb < Eeo for all values of 3D scattering length as- 
According to the Feshbach scheme, a resonance would 
occur when the bound state energy of He coincides with 
the continuum threshold of Hg, leading to the resonance 
condition 

Eeb = Ego = ^^-^ — . (33) 

The bound state energy Eeb of He can be found by pro- 
jecting the total wave function onto the excited Hilbert 
space = Pe"^ , and decomposing into the harmonic os- 
cillators of the excited Hilbert space [H, 

J^^a;)$„^ (2/)$„Jv^z), 
(34) 

where $„ (t) are the dimensionless harmonic oscillators as 
given by Eq. ([8]). Substituting the expansion ([34]) into 
the Schrodingcr equation H'i'e = E'i'e yields 

0= J2 Ec»(^"--^)^"(^'2/'^) 

d 

+2Tras6{r)—r ^ ^ c„P„(a;, y, z) (35) 

where 

Pn{x, y, z) = (V^x)$„„ (2/)$„, iVez) (36) 

Here, n is a shorthand notation for {hx, Uy^riz), and -En = 
{nx + l/2)T]+{ny + l/2) + {nz + l/2)e is the eigenenergies 
of 3D harmonic oscillators. To determine the expansion 
coefficients Cn we project Eq. ([55)1 onto state $„^$„j^$„^ 
with arbitrary n^, riy and 11^, leading to 

_ 27ra,7^7^„(0,0,0)C 
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where C is a constant fixed by normalization of tlie wave 
function. The value of C is related to the limiting behav- 
ior of 'l'e(r) for r approaches zero as 



C = 



9 r ^ ^ 



(38) 



Substituting the solution (|37p for coefficients Cn into 
Eq. (|38|) . we can cancel out the constant C and obtain 
and equation which determines the cigenenergy E, 



1 



(39) 



r-i-O 



where 
*e(f,r) 



E 



^ Vfy?P„(0,0,0)7^n(x,;/,z) ^^^^ 



Notice that the summation in \l/e runs over the excited 
modes with Ux + riy > 0, hence this integral represen- 
tation is valid for £ < min[7y, 1]. i.e., for the eigenstate 
with energy less than Eq + min[77, 1]. Notice that, in the 
limiting case of e — > 0, the zero point energy represents 
the open channel threshold Eo{e ^ 0) = Ego- Thus, we 
can safely use the integral representation (|1T|) for £ < 0, 
and determine the Fcshbach resonance condition. 



To perform the summation over 



and Hz, we 



apply the generating function for the products of Hermite 
polynomials [13] 



n 



^{2axy — a x 



(42) 



and £ = E — Eq denotes the energy shifted by the zero 
point energy Eq = {T] + l + e)/2. It should be emphasized 
that the regularization operator drr and the summation 
in Eqs. ((39)) and (|40)) cannot be interchanged, since the 
series is not uniformly convergent. As we can see later, 
this is related to the divergent behavior of at small r, 
which is regularized by dr{r-). To perform the summa- 
tion in Eq. (|40p we express the denominator in terms of 
the following integral [33} 



Upon inserting Eq. (^1]) into and performing the 

summation according to (j42p . we obtain the following 
integral representation for the wave function 



*e(f,r)=V'^(f,r)+V'^i(f,r), 



(43) 



1 



UxTj + Uy + n^e — £ 



^^^-t[n^V+ny+n,e-£) ^ (41) where 



■qe 



dt 



(27r)3/2 7p ^sinh(?7i) sinh(i) sinh(ei) 



exp 



dt- 



(27r)3/2 7o ^ihni(i) 



: exp 



t £ 



rjx 



2 coth(77t) - y coth(t) 



■ coth(et) 



■ coth(ei) 



(44) 
(45) 



Next we investigate the behavior of the wave function 
for small values of r. Note that the integral of iIj]} is well 
behaved at r = and can be carried out analytically, 
leading to 



small t reads 



^l{£,v) 



(27r)3/2 



dt- 



V(2t) 



^3/2 



(47) 

27rr ^ ' 



^"(^,^ = 0) 



t3/2 



dt- 



10 Vl - 
Vv r(-g/2£) 
27rVir(-£/2e+ 1/2)' 



(46) 



On the other hand, ipe divergent in the limit of r — > 0, 
where the main contribution to the integral comes from 
small arguments t. In the leading order we neglect the 
energy dependence of £, and the expansion of (|44|) for 



By adding the two terms together, we note that the wave 
function "ile diverges in the same way as the Green's func- 
tion in a homogeneous space. This can be understood by 
reminding that at short distances the behavior of the 
wave function is mainly determined by the interaction 
between particles, and has negligible dependence on the 
global trapping potential. Substituting the expression of 
wave function VPg into (p9| . the divergent behavior will be 
regularized by the operator dr{r-), hence make no con- 
tribution to the equation. Therefore, we can subtract 
from '^e the right-hand side of ([T7|) . and obtain a simpler 
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FIG. 3: (Color online) Bound state energy E^h of the closed channel Hamiltonian and bound state energy of the full 
Hamiltonian, as calculated for the quasi- ID limit with e = cj^/ij^y (solid). The results of Ech for small but finite e = 0.1 
are also shown as comparison (dashed). The CIR occurs at the position where E^h coincides with the threshold of the ground 
Hamiltonian E^o, as schematically drawn in the figures (dotted). 



expression determining the bound state energy, 



(48) 



where 



dt 



•qee 



t£/2 



^(1 - e-''*)(l - e-*)(l - e-<^*) 



rjee 



t£/2 



1 



(49) 



Notice that this representation is valid for all values of 
£ < 0. In the limiting case oi £ = 0, the second and the 
third terms remove the divergence for t oo and i -> 0, 
respectively, hence guarantee the convergence of ^(0). 

Finally, we turn to the discussion of quasi- ID regime 
by pushing to the limit of e — >■ 0. This corresponds to 
lowering the axial (z) trapping frequency such that it is 
much weaker than the transverse x-y confinement. In 
this limiting case, the main contribution to the integral 



of T comes from small arguments t, and we perform the 
approximation \/l — e""^* « Vei in the determinator to 
obtain 



dt 



,t£/2 



\/i(l-e-''*)(l 



e-*) 
1 



(50) 



Considering the fact that the CIR takes place when the 
eigenenergy E^b = £ + Eq = EgQ, we reach the equation 
for the CIR condition in physical units as 



(51) 



^(0) 



A numerical evaluation of function J-' shows that this con- 
dition predicts exactly the same CIR position as (|18p , as 
one should expect. In particular, in the cylindrically sym- 
metric case of ?7 = 1 , the integral of J- can be performed 
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analytically, leading to 

F{E) = / dt 



g-t(l-£/2) I 



(52) 



where CH(s,a) = X]^o("' + '^)~* denotes the Hur- 
witz Zeta function. Upon substituting the expression 
above into the CIR condition ([5T|) . and noticing that 
Ch(1/2, 1) = C(l/2), we recover the known result of (l24l) . 

In Fig. 3 we show the bound state energy E^b of the 
closed channel Hamiltonian as a function of 3D scatter- 
ing length as /ay for several different choices of rj. In each 
plot, the location of CIR is indicated by dotted lines as 
the eigenenergy crosses the open channel threshold EgQ. 
We also verify that the quasi-lD regime for two interact- 
ing atoms does not require an infinitesimally small value 
of e, but is already realized for e 0.1 as illustrated in 
Fig. 3 as (red) dashed lines. We have also studied the 
case of e = 0.001, which is close to the realistic experi- 
mental setup [lOl, and note an excellent agreement with 
the quasi- ID energy levels, which on the scale of Fig. 3 
are indistinguishable. 

It is also illuminating to investigate the bound state 
energy of the full Hamiltonian as iJi-oi^(r) = i?t5'(r). 
This can be done by decomposing the wave function into 
the complete set of harmonic oscillators 

*«JV^a;)$„„(2/)$„^(Viz), 



71 .Till ^TI'Z 



(53) 



and following the same procedure as in derivation of (|49| . 
As a result, the bound state energy is determined by the 
equation 



dt 



rjee 



^(1 - e-''*)(l - e-*)(l - e-<^*) 
1 



(54) 



Upon comparing with the equation for the bound state 
energy of the closed channel Hamiltonian ([55]) , we notice 
that in the cylindrically symmetric case of 77 = 1, the 
bound state energy Eeh of the closed channel Hamiltonian 
and the bound state energy Eh of the full Hamiltonian is 
related via 



E, 



eb 



Eh 



for T] = 1. 



(55) 



For general cases of arbitrary 77, the difference between 
Eeh and Eh is in general has dependence of shown 
in Fig. 3 (green/gray solid lines). 



IV. THE TWO-CHANNEL CALCULATION FOR 
GENERAL FORMS OF INTERACTION 

In the previous discussion of CIR, we have employed a 
single-channel zero range pseudo-potential to model the 



3D interatomic interaction. Nevertheless, in the realis- 
tic cases, the atom-atom interactions have finite ranges 
and multi channels with respect to different hyperfine 
states. The 3D scattering length as is determined by the 
inter-channel Feshbach resonance, which is controlled by 
magnetic field B through the _B-dependentinter-channel 
detuning. To take into account of these characters, in 
this section we provide a formal calculation of the effec- 
tive quasi-lD scattering length aiD with a general two- 
channel model which is more realistic for the interatomic 
interaction. Our formal calculation can explicitly give the 
qualitative behavior of oid with respect to the magnetic 
field B, and show clearly that there is only one resonant 
point where aiD takes zero value. This is consistent with 
our above results based on the usage of pseudo-potential. 

In the two-channel model, the atoms are initially local- 
ized in the open channel with a two-atom hyperfine state 
\o)s- During the scattering process, the atomic motion 
in the open channel is assumed to be near resonant with 
a bound state |$) in the closed channel, whose threshold 
energy is much higher than the scattering energy. The 
inter-channel coupling can induce the quantum transition 
between the two channels and then vary the scattering 
amplitude. 

In our calculation, the Hamiltonian for the atomic rel- 
ative motion takes the form 



H,,i^f, + H^+eiB)\c)s{c\ + V. 



(56) 



where e [B) is the threshold energy of the closed chan- 
nel with respect to the hyperfine state |c)s, and can be 
tuned as a linear function of the magnetic field B. In 
this Hamiltonian, is the kinetic energy in the z direc- 
tion, _ff_L is the Hamiltonian of the transverse degrees of 
freedom. As in the above sections, we denote the ground 
state of the transverse Hamiltonian H± as |00)_l. The 
two-channel interaction term V can be further expressed 
as 



V = 14g(r)|o),(o| + yei(r)|c).(c| 
+ [W {r)\o)s{c\ + h.c: 



(57) 



where r is the 3D relative coordinate between the two 
atoms, Vbg (r) and Vc\ [r) are the interatomic interaction 
potentials in the open and closed channels, respectively, 
W (r) the inter-channel coupling, and h.c. denotes Her- 
mitian conjugate. For simplicity, we assume Vbg (f ) is a 
short-range potential with 



Vbg {r) = 0, r > n 



(58) 



Such a simple model is reasonable for the formal cal- 
culation if the realistic atom-atom interaction potential 
is a short-range potential (e.g., the Yukawa potential) 
which decays faster than any power-law function of the 
atom-atom distance r. In the ultracold gases, the atoms 
experience a long-range van der Waals potential which 
behaves as l/r® in the long-range region, rather than a 
short-range one. However, it has been proved that in 
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the low-energy limit there is no difference between the 
behaviors of the 3D s-wave scattering amplitude from 
the van der Waals potential and the short-range poten- 
tial [ssl . [36j . Therefore, we use Eq. (|58p clS cl reasonable 
simplification for the discussion on the low-energy scat- 
tering. 

In our system the spatial parity along the z direction 
is conserved. Since the s-wave bound state in the closed 
channel has an even parity, the Feshbach resonance be- 
tween the open and closed channels can only change the 
scattering amplitude for the states with even parity in the 
open channel. Namely, the scattering amplitude with re- 
spect to odd parity is independent on the magnetic field 
B. Therefore, in this section we only concentrate on the 
scattering amplitude for the states with even parity. 

The ID effective scattering length aiD can be evaluated 
from the analysis of scattering states. Before the scatter- 
ing process, the two atoms are moving in the ground state 
|00) j_ of the transverse Hamiltonian /f_L with relative mo- 
mentum k. Then the relevant scattering state (/c)) 
is given by the Lippmman-Schwinger equation |37| 

|^(+) (fc)) = |fc).|00)^|o). + G^+)(fc2)^|^(+) ^^y^ (59) 

where \\i)z is the eigenstate of with even parity 
z{z\k)z = cosfcz, and g'^^^ (k?) is the free Green's op- 
erator 



^0^ 



T^-Hi-e(B)|c),(c| 



(60) 



Throughout this section, we use the natural unit with 
h = = \. The state |?/;^+^ (fc)) can be rewritten as 



IV'W (fc)) -|^v(fc))|o), + |V'c(fc))k 



(61) 



where the states |^o.c(fc)) are the spatial states corre- 
spond to the open and closed channels, respectively. As 
shown in Appendix \^ the Lippmman-Schwinger equa- 
tion dSSl) of (fc)) can be re-expressed as [38| 

l^o(fc)) = |^[+)(fc))+G(+)(fc)W^|^c(fc)), (62) 
IV'c(fc)) = , \ ,A ^){m\^om, (63) 



/s2 + A (B) ' 



where IV'bg'' C*-)) ^^e background scattering state de- 
fined by the background Lippmman-Schwinger equation 

IV'S^ (fc)) = |fc).|00)^ -f (fc) l^bg (r) |fc).|00)j.. (64) 
Here, g')^^ ik) is the background Green's operator given 

by 



(65) 



fc2 + zO+ - T, + ffi + lig 
In Eq. (|63| the detuning A (_B) is defined as 

A (B) = - £ (B) , (66) 



where E\, is the binding energy of the bound state | $) in 
the closed channel. 

Substituting Eq. into and multiplying the 

Dirac bra on both the left-hand and right-hand 

sides of (|62p . we obtain a linear equation for the factor 
($114^1-00 (fc))- The solution of this equation leads the 
open-channel component of the scattering state 



IV'o(fc)) = lC^ {k)) 



bg 



A(B)-($|iyG[,+^ (k)W\^) 



(67) 



The ID scattering length oid []3) can be extracted 
from the asymptotic behavior of the open-channel spa- 
tial state IV'o (fc)) in the long-distance limit \z\ — > oo. 
As shown in Appendix |D1 in this region we have 

M^it^ik)) = (cos + (A;) e^'-^l^l) 100)^,(68) 



^Git^ (k) 



2k 



where IV'bg'' (k)) is the background scattering state with 
in-going boundary condition and /^^n (k) is the back- 
ground scattering amplitude with low-energy behavior 



fcvcn (k ^ 0) 



-1 -I- ifca^p. 



(70) 



Here, a^fy G real is the background scattering length. 

We can also prove that g[,+' (fc = 0) is a Hermitian 
operator as shown in Appendix [d1 Therefore, in the 
long-distance region we have 

,(z|Vi+Ufc)> = (cosfcz + /even (fc;S)e^'l"l) |00)i (71) 

with the scattering amplitude given by 

./cvcn(fc^O;B) = -l + ifcaiD(B) (72) 

in the low-energy limit. Here, the _B-dependent scatter- 
ing length fliD (B) takes the form 



fliD (B) = a^g + lim 



1 



[m\4V (k)) 



A{B) 



(73) 



To obtain the expression (1751) . we have also used the re- 
suits \ijll\k)) (xk and \^pi-\k)) = which 
are proved in the limiting case of fc — > as shown in 
Appendix |DJ 

The result of aiuiB) expressed in (|73p is the central 
result of this section. It gives the behavior of oid as a 
function of the magnetic field B through the magnetic 
Feshbach resonance. Notice that A{B) is a linear func- 
tion of B, thus one can immediately read from Eq. (j73p 
that there is only one resonant point 

B. = ^ {e, - ($|M/G+ (0)1^1$) + AJa%) (74) 
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leading to zero value of cid (B*) = 0. Here we have 
assumed e{B) = fiB and is defined as 



(75) 



This observation is in consistent with our findings in the 
previous sections by assuming a S pseudo-potential. We 
stress that this formal scattering calculation is based on 
a very general form of the quasi- ID Hamiltonian (|56p . 
and the qualitative conclusion of one single CIR remains 
valid disregard to the details of transverse confinement 
or the interatomic potential. 



V. CONCLUSION 

In summary, we have presented a detailed analysis of 
the scattering process of two interacting atoms confined 
in a quasi-lD harmonic trap. We analyzed the possibility 
of confinement- induced resonance (CIR), where the ID 
scattering degenerates to a total reficction as a ID gas of 
impenetrable bosons. By modeling the interatomic inter- 
action by a 3D S pseudo-potential, we have shown that 
only one single CIR would be present in a quasi-lD trap 
with arbitrary transverse anisotropy. The location of this 
resonance has non-monotonic dependence on the trans- 
verse anisotropy ratio t] = ujx/ujy. For slight transverse 
anisotropy with 77 ^ 1, the position of CIR varies with the 
same trend as the experimental observation of CIR2 [s^] , 
but deviates from the measurements in a quantitative 
level. The other resonance observed in experiment, CIRl, 
cannot be accounted for the present theory. To extend 
the discussion to general forms of s-wave interaction, we 
presented a formal theory within a two-channel model, 
and draw the same qualitative conclusion of one CIR for 
transverse confinement with arbitrary anisotropy. 

After finalizing this manuscript, we are aware of a work 
recently posted by Peng et al. [39| , which discusses CIRs 
in quasi- ID geometries with transverse anisotropy and 
reaches similar results as our calculation with pseudo- 
potential. However, the formal theory which is valid for 
a more general form of interaction is not discussed in |39( . 



Appendix A: Background Lippmman-Schwinger 
Equation 

In this appendix we derive the Lippmman-Schwinger 
equations (p^ and for the scattering state li/i^"*") (fc)). 
We first substitute Eq. ((6T|) into ((59)) . and get 



1^0 (fc)) = |fc).|oo)i + G[+'(fc2)ybg|Vo(fc)) 



GL+'(fc')M^|Vc(A:)) 



(Al) 



\^c{k)) = Gi,{e)V,x\iP,{k))+Gi,{e)W\i^o{k)) 

(A2) 

with the free Green's operators G^^^^ and Gfc in the open 
can closed channels, respectively. 



G'^\k') = 

Gfe(fc') = 



1 



fc2 + iQ+ -T^- Hj_ 
1 



fc2 -T^-Hj_- e{B) 
By defining the closed channel Green's function Gci 

Gci(fc') 



(A3) 
(A4) 



(A5) 



k^-T,-H^~ Vci - e(S) 
we obtain the relations 

G\:\k') = GW(fc^)-GW(fc^)FO^^)GW(fc^); 

(A6) 

Gfe(fc') = Gei(fc2)-Gci(fc2)V^(^')Gf,(fc2), (A7) 
where G'^\k'^) is defined in Substituting Eqs. 



[6]) and (jA7[) into the last terms of the right-hand-side 
(r.h.s.) of Eqs. (jAip and ()A2|) . and using the Lippmman- 
Schwinger equation (|64p for the background scattering 

state l^bg' (^)) ^^^^ approximation 



Gci(fc' 



|$)(cf>| 



k^ + A {B) 



(A8) 



with A(i?) defined in Eq. we get equations (p^ 

and (|63p for the states \il)o [k]) and \ipc {k))- 



Appendix B: Low-energy Behavior of ID Scattering 
State and Green's Operator 
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In this appendix we discuss the low-energy behavior 
of the scattering state in a single-channel ID scattering 
problem. The results of this discussion are not directly 
used in the main text of our paper, but they are prereq- 
uisite for the calculation in Appendix [Dl 

We consider a ID scattering problem with symmetric 
short-range potential with 



U{z) 
U{\z\ > r,) 



0. 



(Bl) 
(B2) 
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The wave function (fc; z) of the scattering state 
with even parity hence satisfies the Lippmman-Schwinger 
equation 



0(+)(fc;z) = cosfcz+ / dz' (fc; z, z') 



(B3) 



where ^q"'''' (fc; z, z') is the ID Green's function 



ffo^' {k;z,z') = ^{z\ 



1 



2k 



(B4) 



It is straightforward to show that the Lippmman- 
Schwinger equation (jB3[) is equivalent with the time- 
independent Schrodingcr equation [37| 



dz2 



(fc; z) + U (z) (A:; z) = (y^. 2)(B5) 



with the boundary condition 

(/c; \z\ > r,) cosfcz + /(A:)e*'=l^l. (B6) 

Here, f{k) is the scattering amphtude. 

To investigate the low-energy behavior of (p^^^ (k; z) 
and / (fc), we expand t/)'^^ (fc; z) in series of k 

</)W(fc;z) = (j)^{z) + (t)ir{z)k + i(t)u{z)k + 0{k'^) , 

(B7) 

where (z) and (pn (z) are assumed to be real functions 
of z. Setting fc = in Eq. (jBSp . we get the equation for 

4>o {z) 



-^(j^o {z) + U{z)cj)o (z)-0. 



(B8) 



Considering the fact that the potential U (z) takes zero 
value in the region of |z| > r*, the solution (f>o{z) of the 
equation above would be a constant within this outer 
region. Therefore, in the inner region of |z| < r*, we 
acquire the following boundary condition 



dz 



00 (z) 



0. 



(B9) 



\z\=r. 



This boundary condition is a strong requirement such 
that only some special kinds of U (z) (e.g., the most triv- 
ial case with U (z) = 0) can lead to a non-zero solution 
of Eq. (jBSp . Here we assume our potential are not of 
this special form, then we have 



00 (z) = 0, 
(fc 0;z) cx k. 



(BIO) 



Substituting the equation above into (jB6[) . we obtain the 
low-energy behavior 



/(fc->0) = -l + C'(fc) 



(Bll) 



of the scattering amplitude with even parity. 

Now we focus on the linear terms of k. Expanding 
Eqs. (jBSP to the first order of k, we find that the func- 
tions 0ir (z) and (z) defined in Eq. ()B7| also satisfy 
Eq. (jB8[) . The expansion of Eq. (|B6p further shows 
that (z) takes constant value in the outer region of 
|z| > r*. Thus, in the inner region of |z| < r* the wave 
function (z) satisfies the same equation and boundary 
condition as (z), leading to 



01r (z) = 00 (z) = 0. 



(B12) 



As a consequence, the scattering state 4>{k] z) with even 
parity acquires the following form in the low-energy limit 
fc -> 

{k 0; z) = ik(f)u (z) ; 0h (z) e real (B13) 
which yields 

f{k^0) = -l + ika + O{k^); a G real. (B14) 

Next, we discuss the low-energy behavior of the ID 
Green's operator 



1 



k^-T,-U {z) 



(B15) 



In the limiting case of fc — > 0, the Green's operator can 
be expressed as 



,(z|g(+)(fc2^0)|z') 



dp 



c^(+>{p;z)c^(+Hp;z') 



p^ 



The convergence of this integral representation is guaran- 
teed by noticing that 0"*" {p; 2) oc p in the limit of p — 0, 
as wc have shown above. Besides, the expression above 
also satisfies 

,(z|gW(fc2 -> 0)|z% = .(z'|g(+)(A:2 ^ 0)|z):, (B16) 

which indicates that the Green's operator converges to a 
Hermitian operator in this low-energy limit. 

Finally, we discuss the low-energy behavior of the scat- 
tering state (fc, z) with in-going boundary condi- 
tions. The Lippmman-Schwinger equation for 0^"' (fc, z) 
reads 



(j)'- •* (fc, z) = cos kz ■ 



dz' 



'gi-hk;z,z') 



xC/(z')0(^Ufc,z')l, (B17) 



where (fc; z, z') is defined as 

5o~'' {k]z,z') = 



fc2 - iO+ - 



2k 



(B18) 
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It is easy to prove that, 0^") (fc, z) and (j)^^^ (fc, z) sat- 
isfy the same ID Schrodinger equation with the same 
eigenvalue, and they also share the same parity. Thus, 
(f>^~^ (fc, z) is proportional to t/)*^"'"-' (fc, z) and we have 



(fc, z) = ^ (fc) 



By expanding Eqs. ()B3P and ()B17P to the zeroth order 
of fc, and using Eqs. (HH, ((BT3)) . (|BT8)) and ([STOt . we 
obtain the following equations 



I- -J dz'^u {z') U {z') = 0; (B20) 
l + ^e(fc->0) f dz'(l^u{z')Uiz') = 0, (B21) 



which leads to 



e(fc^o) = -i 



(fc,z) = (fc,z) 

in the limit of fc ^ 0. 



(B22) 
(B23) 



Appendix C: The Asymptotic Behavior of ID 
Green's Operator 

In this appendix we discuss the asymptotic behavior of 
the ID Green's operator. We focus on the case with even 
parity since it is the only relevant term in our problem. 
In this case, the ID Green's operator reads 



^Gvcn (^i Z, Z ) 



P +iO+ -T^-U (z) 



-Povonl^ )2,(C1) 



where Pcvon is the projection operator to the space of 
states with even parity of z. 

We first investigate the asymptotic behavior of the free 
even Green's function 



5cvcn0 Z, Z ) — -PcvGn<7o ^ (^i Z , Z ) 



(+) 



(C2) 



with ' {k; z, z') defined in (jB4[) . In the case of \z\ > 
and jz'l < r*, we can show clearly that 



,Q (fc; z, z') = -i^e'''l^l cos kz'. 



(C3) 



Now wc consider the behavior oi gitcn (k; z, z'\ We no- 
tice that (7oven (fc; z, z') satisfies the Lippmman-Schwinger 
equation 



+ / d^"5*:in0 (fc; ^") V iz") g^+>^ (fc; z" , z') (C4) 



Therefore, according to Eq. (jCSp . we know that in the 
region with \z\ > and \z'\ < we have 



g(+)„(fc;z,z')=-*^e^'l^lx(fc,z')*, (C5) 



(B19) where 



x{k,z')* = coskz' + J dz" cos {kz") 

xUiz")9^otLik;z'\z')]. (C6) 

Here, we have used the fact that U (z) = in the region 
\z\ > r^. Notice that Eq. (jC6p implies 



X{k,z') = cosfcz' + J dz^U{z") 

xgU,Ak;z",z') cos kz"], (C7) 
where the Green's function given {k; z, z') is defined as 

fcvcn (^i Z, Z ) 



P -iO+ -T^-U (z) 



3%.(C8) 



It is apparent that the right-hand-side of Eq. (|C7[) is 
equivalent with the Lippmman-Schwinger equation (|B17[) 
for the scattering state (j)~{k,z) with in-going boundary 
condition, leading to xik-,z) = 0~(fc, z). Thus, in the 
region |z| > wc have 

gitL (fc; z, z') = -z^e''=l-l0(-) (fc, z'r ■ (C9) 



Appendix D: Background Scattering State 

In this appendix we investigate the low energy behav- 
ior of the background scattering state IV'bg^ ik))\o)s, the 
background Green's function g[^^ and the background 



scattering amplitude Z^^^. The state I'/'bg {k)) is given 
by the Lippmman-Schwinger equation (j64p . For conve- 
nience we define the states l^'bg'' (k)) p.Q a-s 



\4t^ik))p - (fc)), 



(Dl) 



= (D2) 

where = |00)_l(00| is the projection operator to the 
transversally ground state. Following a similar procedure 
as in Appendix|Xl we can write the Lippmman-Schwinger 
equation ((M)) as 



\^iVik))p = |0<.+'(fc)).|OO)^ 



~G 



(+) 



p (fc)xiv[,:^(fc))c 



iV'i? (fc)>Q = G[3+)(fc)X|^(+)(fc))j 



(D3) 
(D4) 
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where 



G^^^ (fc) = 



1 



1 



(D5) 



P -T,+ [l- Pg) [H^ + Fbg j (1 - 

(D6) 

1 = (h^ + Kg) (l - Pg) + h.c, (D7) 

and |(/'^'' is the ID scattering state with respect to 
the potential _L(00|Vbg (t") |00)x in the z direction. We 
further rewrite Eq. (|D4p as 

(D8) 



with 



d,: = 



fc2 



(D9) 



Here, |ci) is the eigenstate of the Hamihonian + 
(^1 - Pg) (^Hi_ + Vbg) (^1 - Pg) with eigenenergy e^. 

Combining Eqs. (|D3IID4p and (|D8IID9p together, we get 
a set of equations for the parameters di, leading to the 
following expression for the state \ipl'^^ {k))p 

+ ^ G(.+) (fc) X|c,)i?-.i (fc) (c,|l|0(.+) (fc)).|00)^, 



(DIO) 



where i^^^ (fc) are the elements of the inverse of matrix 
Dik) 

D,, (k) = (fc2 _ £,.) 5,.^. _ (q|XG^+) (fc) ^|c,>. (Dll) 

Now we investigate the behavior of IV'bg' (k)) p in the 
limit of fc — > 0. As shown in Appendix iBl in this low en- 
ergy limit the ID scattering state \4>'p '' (k))z is a linear 

function of k, and the Green's function G^'' (fc) becomes 
a constant Hermitian operator. Substituting these prop- 
erties into Eq. (jDlOp . we get 



IV-i? {k ^ 0))p cx k. 



(D12) 



That is, the background scattering state IV'bg^ {k))p also 
linearly depends on k in the low energy limit. Substitut- 
ing (|D12p into (|D4p . we get the low energy behavior of 
the total background scattering state 



\^iV {k ^ 0)) (X k. 



(D13) 



Using this result, and following the same analysis as in 
Appendix [B] wc can show that the background Green's 



operator g\^^ (k) defined in ([65)) converges to a Hermi- 
tian operator in the low-energy limit fc — >■ 0. 

In the above we obtained the low-energy behavior 
of the background scattering state and the background 
Green's function. Next we consider the low-energy be- 
havior of the background scattering amplitude. We also 
assume Vbg (r) = when r > r^. Then in the region 
\z\ — )• oo, the wave function of IV'bg^ decays to zero 

while the wave function of I^Atg^ ik))p can be expanded 

as 

±(OOU(z|vS^(fc))p =cosfcz + /X(fc)e''l^l (D14) 
with the background scattering amplitude 



fcvon (k) — fcvcn (k) * 



1 



(0(r^(fc)U(OO|f/|c,) 



2k ^ ^^'^^ 



xD-Mk) {c,\U\^^+'> (fc)),|00). 



(D15) 



Here, f^^n (k) is the scattering amplitude for the state 
\4'p'^ {k))z, and 10^^ (fc))z is the ID scattering state with 
respect to j_(00|Vi3g (r) |00)_l under the in- going bound- 
ary condition. To derive the equation above, we have 
used the behavior of the ID Green's function G^"* (k) in 
the region z > (see Appendix [C|) 

,(2U(00|G(,+)(fc)|00)i|z% = 

-*^e^'=l^l.(0(r^(fc)|/).. (D16) 



As we have shown in Appendix [B] the scattering am- 

4" 



phtude /^on {k) and scattering states \4)"p^ {k))z have the 



following properties in the limit fc — > 

— 1 + ika^Y), afn G real; (D17) 



/even {k) 



-I0P ^ (fc)). 



(D18) 



Then it is obvious to sec that Eq. (jD15|) can be reformed 
as 



fcvcn (k) — 1 + ika^f), 
where the background scattering length 



(D19) 



^ID 



*1D 



xDT-} (fc) {c,\U\cj)''+^ (fc)>.|00)i (D20) 

takes real value. 

In the end of this appendix, we stress that the analysis 
in Appendix [C] on the asymptotic behavior of the ID 
Green's operator can be straightforwardly generalized to 
the quasi- ID case and we have 

(fc) = -^^mM~J (k) I (D21) 
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in the region of |z| > r*. Here, the scattering sate with 
l^bg^ (fc)) with in-going boundary condition is given by 



1 



Gi~J (k) = ^ ^ . (D23) 

Ht,Hk)) = \k).m^ + Gl-Hk)V^^{r)\k),m^ k^-iO+~T,+H^ + V^, 

(D22) 
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